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The paramagnetic metallic phase of the one-dimensional 
Kondo lattice model is studied by the density matrix renor- 
malization group method. We observe charge and spin Friedel 
oscillations. They reflect the long range charge-charge and 
spin-spin correlation functions. The observed oscillations are 
consistent with a Tomonaga-Luttinger liquid. From the pe- 
riod of the oscillations it is concluded that the Fermi surface 
is large, including both the conduction electrons and the lo- 
calized spins, kF = 7r(l + n c )/2, where n c is the density of 
conduction electrons. 

PACS numbers: 75.30. Mb, 71.27.+a, 71.28.+d 



Electron density oscillations as a response to a local 
perturbation provide a clue to understanding electronic 
states of a system. Response to an impurity potential 
is known as Friedel oscillations. The period of the os- 
cillations is given by the diameter of the Fermi surface 
which is defined by the singular points in the momentum 
distribution function. As a local perturbation it is also 
possible to use a magnetic impurity. Then spin density 
oscillations are induced. Since these oscillations are the 
origin of the Rudcrman-Kittel-Kasuya-Yosida (RKKY) 
interaction, the Friedel oscillations are responsible for the 
magnetically ordered structure of rare earth metals. 

Generally speaking, a key to understanding the physics 
of heavy fermion systems is the competition between the 
RKKY interaction and the Kondo screening effect of the 
magnetic moments. The Kondo lattice model is one of 
the canonical models for heavy fermions and much ef- 
fort has been devoted to elucidate the properties of the 
model. Recently the ground state phase diagram of the 
Kondo lattice model was completed in one-dimension 
[0,1). There are three phases in one dimension: a fer- 
romagnetic metallic, a paramagnetic metallic, and an in- 
sulating spin liquid phase. 

Our concern in the present paper is the metallic phase 
of the one-dimensional Kondo lattice model. In one 
dimension it has been established that most interact- 
ing metallic systems belong to the universality class of 



Tomonaga-Luttinger liquids ||. The asymptotic forms 
of charge- and spin-correlation functions are 



(n(x)n(0)) = K p /(nx) 2 + A x cos(2k F x)x~ 1 ~ 
+A 2 cos(4:k F x)x~ 4:K >' , 
(S(x) ■ S(0)) = l/(irx) 2 + B 1 cos(2k F x)x- 1 - K " 



(1) 

(2) 



where hp — irp/2, with p being the density of charge 
carriers, is the Fermi momentum and K p is the corre- 
lation exponent B. In the above equations logarithmic 
corrections to the 2kp correlations have been omitted. 
The momentum distribution function around kp shows a 
power law singularity ~ 1/2 — sgn{k — kp)\k — kp\ a 
with a = {K p + X/Kp — 2)/4. The anomalous power law 
decays of the correlation functions naturally reflect them- 
selves in the Friedel oscillations: the asymptotic form of 
the charge density oscillations induced by an impurity 
potential is 



8p(x) ~ Ci cos(2kpx)x 



(-l-K p )/2 



Ci cos(Akpx)x 



-2K„ 



(3) 



as a function of the distance x from the impurity 
Analogously, the spin density oscillations induced by a 
local magnetic field behave as 



<j(x) ~ D\ cos{2kpx)x' 



-K a 



(4) 



Concerning the one-dimensional Kondo lattice model 
the paramagnetic metallic state is expected to belong 
to the class of Luttinger liquids. However in this case 
the position of the Fermi points is already a nontrivial 
problem since the model consists of two components with 
completely different characters, the conduction electrons 
and the localized spins. 

There are two different points of view concerning the 
above question. If the interaction between the conduc- 
tion electrons and the localized spins is strictly zero, it 
is clear that the singularity in the momentum distribu- 
tion function of the conduction electrons is determined 
only by the number of conduction electrons. Thus, if 
the singular points are not affected by the interaction, 
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kp = im c /2 is expected. On the other hand, a differ- 
ent answer is obtained by identifying the Kondo lattice 
model as an effective model for the periodic Anderson 
model. In the periodic Anderson model the conduction 
electrons and the / electrons are mixed with each other 
through the hybridization matrix elements. Therefore it 
is naturally expected that Uf is determined by the to- 
tal density of both the conduction electrons and the / 
electrons; kp — 7r(l + n c )/2. 

In order to draw a conclusion on whether the Luttinger 
sum rule includes localized spins or not, it is necessary 
to deal with large-size systems because the Luttinger liq- 
uid is characterized by the long range correlations and 
the singularity in the momentum distribution function is 
clearly defined only in the infinite system. For this pur- 
pose the density matrix renormalization group (DMRG) 
method developed by White [0 is most promising since in 
the DMRG we can study long chains and obtain results 
with only small systematic errors, which can be estimated 
from the eigenvalues of the density matrix. 

In this paper we calculate the Friedel oscillations of 
the one-dimensional Kondo lattice model by using the 
DMRG. The charge density oscillations are induced nat- 
urally by open boundary conditions and the spin den- 
sity oscillations are introduced by applying local mag- 
netic fields at both ends. From the period of the Friedel 
oscillations we can determine the Fermi momentum. Sur- 
prisingly the Friedel oscillations are compatible with the 
large Fermi surface kp = tt(1 + n c )/2 in spite of the fact 
that the charge degrees of freedom are completely sup- 
pressed for the / electrons in the Kondo lattice model. 

The Hamiltonian we use in the present study is the 
usual one-dimensional Kondo lattice model, 

H = -*5>*Wi* + H.c.) + jJ2 S i s i ( 5 ) 

where cl (ci a ) is the creation (annihilation) operator 
of a conduction electron at the i-th. site, and s' 1 = 

( l l 2 ) Y.aa< 4a T aai c ™'i with T ^ = x,y,z) being the 
Pauli matrices, are the spin density operators of the con- 
duction electrons. The spin densities are coupled to the 
localized spins through an antiferromagnctic exchange 
coupling J. 

Here we briefly summarize the phase diagram of this 
model At half filling, n c = 1, the ground state is a 

spin singlet with a gap for excitations [P-|l2t . Away from 
half filling the ground state is ferromagnetic in the strong 
coupling limit j|13| , pT[ . The ferromagnetic ground state is 
continuously connected to the low carrier density limit, 
n c — > 0, where the ferromagnetic state survives even in 
the limit of weak exchange couplings. The paramagnetic 
metallic ground state exists only for weak couplings for 
finite carrier densities. The phase boundary between the 
ferromagnetic state and the paramagnetic state had been 
determined by the numerical exact diagonalization of fi- 
nite clusters, and it was shown that the critical value J c 
goes up with increasing the carrier density ||. 



To treat large systems with sufficient accuracy, small 
truncation errors in the DMRG calculation are necessary. 
Since the truncation errors decrease with increasing J, we 
chose n c = 4/5 for the calculation of Friedel oscillations, 
where the paramagnetic state remains up to J c = 3M. 
At this filling there is a small additional ferromagnetic 
region below J c — 3.0t. It has already been reported in 
previous exact diagonalization studies at slightly different 
carrier concentration n c = 0.75 Our DMRG results 
for clusters of size N = 10 and 20 confirm the existence of 
this ferromagnetic region between J = 1.6t and J = 1.8t 
at n c — 4/5. The total spin of the ferromagnetic state 
is the same as that for strong coupling, and the phase 
transitions are caused by simple level crossings between 
the two lowest states with S = and S = N(l — n c )/2. 
We calculate Friedel oscillations in the two paramagnetic 
states at J = 2.5t and J = 1.5t. 

We first consider the paramagnetic state at J = 2.5t. 
The charge density oscillations induced by the open 
boundary conditions are presented in Fig. 1 by the solid 
line. As expected, long range oscillations characteristic 
of a Luttinger liquid are induced by the boundary con- 
ditions. Since the weak decay of the oscillations makes 
it difficult to determine the correlation exponents K p by 
the present system size, we focus our attention mainly on 
the period of the oscillations. The Fourier components 
of the oscillations, presented in Fig. 2, show a clear sin- 
gle peak at q = 2tt/5. The single peak at this wave 
number is natural in the strong coupling limit where the 
conduction electrons and the localized spins form local 
singlets, leading to complete spin-charge separation. In 
this case we can treat the charge part as spinless fermions 
whose correlation functions are characterized by its sin- 
gle 2kp structure given by 27r(l — n c ). The peak at 
q = 27r/5 corresponds to the 2kp structure of the spin- 
less fermions which is equivalent to 4/cf structure of the 
original fermions. The result of Fig. 1 means that the 
charge density oscillations at J = 2.5t are already well 
characterized by the nature of the strong coupling region. 

From the above results, it is not possible to draw a con- 
clusion about whether the Fermi surface includes the lo- 
calized spins, because 4kp — 27r(l+n c ) = 2irn c (mod2n). 
In order to see the contribution of the localized spins we 
next calculate spin density oscillations. If the 2kp struc- 
ture is found in the spin density Friedel oscillations re- 
flecting the spin-spin correlations of Luttinger liquids, it 
may be possible to argue whether the Fermi surface in- 
cludes the localized spins or not: if the localized spins 
contribute to the Fermi surface a 2kp — 7r/5 structure 
should appear but if it does not a 47r/5 structure appears. 
To induce the Friedel oscillations of spin density, we ap- 
ply local magnetic fields, Hi oca i — 2h(S* — sf — + s^) 
to the spins and conduction electrons at the boundary 
sites. These local magnetic fields with opposite direc- 
tions for the two boundary spins induce oscillations that 
are odd with respect to reflection. This feature has been 
used to confirm the convergence of the DMRG. 

The spin density oscillations are presented in Fig. 3 by 
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the solid line and their Fourier components are shown 
in Fig. 4. We clearly see long range oscillations with 
q = 7r/5. They are consistent with a Luttinger liquid 
prediction, Eq. (4), if we assume the large Fermi surface 
where both the conduction electrons and the localized 
spins contribute to the Fermi volume: 



n c ). 



(6) 



Now we proceed to the paramagnetic metallic phase 
in the weak coupling region. The result of the charge 
density oscillations at J = l.bt is presented in Fig. 1 by 
the broken line. In this case too long range oscillations 
are induced and their period is the same as J = 2.5t. The 
similarity is also seen in the Fourier components shown 
in Fig. 2 by the broken line. For these coupling strengths 
we cannot see a peak at q = 2kp = tt/5. This means the 
amplitude of the 2k p oscillations of the charge density is 
still negligible even in the case of J = 1.5t. However, as is 
shown in the inset of Fig. 2 we find clear 2kp oscillations 
in addition to the dominant Akp oscillations at a smaller 
coupling J = l.Ot which is consistent with the general 
form of the Luttinger liquids. 

For the spin density at J = 1.5t, the oscillations and 
their Fourier components are presented in Fig. 3 and Fig. 
4, respectively, by the broken lines. Although a small 
structure is also found at q — 37r/5, the dominant com- 
ponent is at q = tt/5. The structure at q — 3tt/5 is 
considered to be induced by the coupling mode of the 
spin density oscillations with q = tt/5 and the charge 
density oscillations with q = 2n/5 whose amplitude is 
larger than that of J = 2.5t. Since the structure at 
q = 37r/5 decreases rapidly with an increasing system 
size compared with the structure at q — tt/5, we observe 
that it is not an intrinsic property of the infinite system. 
Combined with the results of the charge density oscilla- 
tions it is concluded that the paramagnetic state at 1.5f is 
also characterized by the Luttinger liquid with the large 
Fermi surface. 

At representative values for the exchange coupling con- 
stant J = 2.5/j and J = 1.5t, we have seen that the 
paramagnetic metallic phases of the Kondo lattice model 
show Friedel oscillations characteristic of a Luttinger liq- 
uid with a large Fermi surface. Similar results are ob- 
tained at different concentrations of conduction electrons, 
n c = 2/3, J = 2.0/j and n c = 6/7, J = 1.7t. These re- 
sults are consistent with the previous work on the t-t' 
Kondo lattice model for which it was shown exactly that 
its strong coupling limit is described by a Luttinger liq- 
uid with the large Fermi surface jll) . Therefore it is nat- 
ural to conclude that the paramagnetic metallic phase 
of the Kondo lattice model has a large Fermi surface in 
general. This conclusion is consistent with a variational 
Monte Carlo study using Gutzwiller projected hybridiza- 
tion form jl6) . This conclusion is also consistent with the 
bosonization study by Fujimoto and Kawakami p7| . Al- 
though their work was recently criticized by White and 
Affleck Jia], it should be noted that the model used by 



the latter authors is different from the usual Kondo lat- 
tice model. 

The DMRG method has made it possible for the first 
time to observe the Friedel oscillations corresponding to 
the large Fermi surface in the paramagnetic phase of the 
weak coupling region. However, it is still difficult to see 
them in the weak coupling limit J < t, where the length 
scale of the Luttinger liquid properties clearly exceeds the 
system size studied in the present work. Calculations on 
longer systems and the determination of the correlation 
exponent K p are important areas for future investigation. 
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FIG. 1. Charge density oscillations of the Kondo lattice 
model. The system size is 60 sites and the carrier density 
is ric — 4/5. The solid line and the broken line correspond 
to J = 2.5* and J — 1.5*, respectively. Typical truncation 
errors in the DMRG calculations are 1 x 1(T 6 for J = 2.5* 
and 3 x 1(T 6 for J = 1.5*. 



FIG. 2. Fourier components of the charge density oscilla- 
tions. The system size is 60 sites and the carrier density is 
n c — 4/5. The solid line and the broken line correspond to 
J = 2.5* and J = 1.5*, respectively. The Fourier transfor- 
mation is carried out by using the central 50 sites. The inset 
shows the result for J = 1.0*, for which typical truncation 
errors are 4 x 10~ 5 . 



FIG. 3. Spin density oscillations of the Kondo lattice 
model. The system size is 60 sites and the carrier density 
is n c = 4/5. The solid line and the broken line correspond 
to J = 2.5* and J = 1.5*, respectively. The strength of the 
local magnetic field h is 0.1*. Typical truncation errors in the 
DMRG calculations are 1 x 10~ 6 for J = 2.5* and 3 x 10~ 6 
for J = 1.5*. 

FIG. 4. Fourier components of the spin density oscilla- 
tions. The system size is 60 sites and the carrier density is 
n c = 4/5. The solid line and the broken line correspond to 
J — 2.5* and J = 1.5*, respectively. The Fourier transforma- 
tion is carried out by using the central 50 sites. 
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